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Abstract 

The process pp — > A^A~ is investigated within the handbag approach. It is 
shown that the dominant dynamical mechanism, characterized by the partonic sub- 
process uu — > cc, factorizes in the sense that only the subprocess contains highly 
virtual partons, a gluon to lowest order of perturbative QCD, while the hadronic 
matrix elements embody only soft scales and can be parameterized in terms of 
helicity flip and non-flip generalized parton distributions. Modelling these parton 
distributions by overlaps of light-cone wave functions for the involved baryons we 
are able to predict cross sections and spin correlation parameters for the process of 
interest. 

1 Introduction 

The FAIR project at GSI with the HESR antiproton program will offer ideal possibilities 
to study exclusive channels in pp annihilation. The energy of the antiproton beam suffices 
to produce pairs of heavy baryons, as for instance A+A~, in proton-antiproton collisions. 
Prerequisite for the measurement of such processes is, however, that their cross sections are 
sufficiently large. The purpose of this work is to estimate the cross section for the reaction 
pp — >• A+A^. This estimate is based on a generalization of the, so-called, "handbag 
approach", which has been developed for deeply virtual electroproduction of photons 
and mesons [1]. It has been firmly shown that the amplitudes for these processes are 
described by convolutions of hard subprocess kernels and, so-called, "generalized parton 
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distributions" (GPDs), which encode the soft, non-perturbative physics. A second class 
of hard exclusive reactions to which the handbag approach can be applied to, is formed by 
wide-angle reactions. The hard scale for these reactions, necessary to justify factorization 
into hard and soft physics, is provided by the momentum transfer from the incoming to 
the outgoing hadron or, in other words, by large Mandelstam variables —t and —u, instead 
of the photon virtuality. Although, in contrast to deeply virtual processes, rigorous proofs 
do not exist, arguments for the dominance of the handbag contribution have been given 
in some cases [2, 3, 4]. 

Along the lines of argumentation for widc-angic Compton scattering [2] we demon- 
strate that, under the assumption of restricted parton virtualitics and transverse mo- 
menta, the pp — > A+A^ amplitudes factorize into hard subprocesses and moments of the 
u ^ c transition GPDs. Owing to the relatively heavy charm quarks (with mass rric), 
the process possesses the large intrinsic scale 4m^ ~ 6.3 GeV^ to which the momentum 
transfer adds. The large intrinsic scale leads to the highly welcome consequence that large 
—t and —u are not required for wide-angle Compton scattering. Therefore the 

arguments for factorization of the pp — > A+A~ amplitude already hold for small scattering 
angles. 

It is to be stressed that the process of interest has already been investigated long 
time ago [5] in an approach that bears resemblance to the handbag one. As compared 
to [5] the progress achieved here manifests itself in the more solid theoretical basis. The 
normalization of the cross section is now determined for given GPDs and not fixed by 
symmetry arguments as in [5]. Moreover, the complicated structure of the proton which 
is described in terms of eight GPDs, can be easily taken into account in the handbag 
approach and therefore the spin dependence of the process pp — > A+Aj can readily be 
treated. 

The plan of the paper is the following: The arguments for factorization and the prin- 
ciple structure of the handbag contribution are presented in Sect. 2. In Sect. 3 the ampli- 
tudes of the process pp — > A+Aj are discussed while the subprocess amplitudes are given 
in Sect. 3.1. Various models for the GPDs are presented in Sect. 4 and predictions for 
observables are discussed in Sect. 5. The summary is given in Sect. 6. In two appendices 
various kinematical formulas and the properties of the transition GPDs are compiled. 

2 Factorization 

2.1 The basic idea 

We assume that the graph shown in Fig. 1 dominates the process of interest for Mandel- 
stam s well above the kinematical threshold and in the forward scattering hemisphere. 
We are going to show in this section that the baryons emit and re-absorb soft partons, 
i.e. quarks that are almost on-shell and move almost coUinear with their parent baryon. 
In other words, the contribution from the graph shown in Fig. 1 factorizes into a hard 
partonic subprocess and soft hadronic matrix elements. The arguments for factorization 
go along the same lines as for wide-angle Compton scattering [2], except that here we 
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Figure 1: The handbag contribution to the process pp A+A^ . The momenta and 
hehcities of the baryons and quarks as well as the colors of the latter, are also specified. 

have to consider the double handbag and we have to take care of the unequal mass kine- 
matics For the p ^ transition a w-quark is to be emitted and a c-quark reabsorbed; 
the valence quark content is to be changed. The emission of a light quark other than the 
u one from the proton does not lead to a Ac baryon in the final state. Intrinsic charm 
would allow for subprocesses like 

cu uc. (1) 

This is a higher Fock-state contribution which is expected to be suppressed by general 
arguments. Moreover, according to present PDF analyses, e.g. [7, 8], the charm content 
of the proton is tiny. Even at a scale of, say, 25 GeV^ it only amounts to about 10^^ of 
the M-quark content at x ~ 0.6. Thus, we merely have to take into account the subprocess 

uu cc. (2) 

This subprocess goes along with p —>■ A+ transition matrix elements which, as we are 
going to show below, can be parameterized by m — > c GPDs. One may also think of time- 
like matrix elements or GPDs representing pp —>■ uu and cc — > A+A~ transitions. Such 
time-like GPDs have been employed in the analysis of wide-angle time-like processes as, 
for instance, two-photon annihilation into a pair of hadrons [9, 10, 11]. To lowest order of 
QCD the subprocess is mediated by one-gluon exchange which evidently forbids the time- 
like reaction mechanism by color conservation to this order of accuracy. It is, however, 
allowed to higher orders of QCD although, beyond representing an correction, expected 
to be smaller than the mechanism depicted in Fig. 1 since the p A^ transitions seem 
to favored over the time-like ones. We therefore do not consider the time-like mechanism. 

A u c transition GPD is a function of three variables, a momentum fraction x 
defined by the ratio of light-cone plus components of the average parton and baryon 

^With the exception of the latter issue this is similar to proton-proton elastic scattering which has 
been briefly discussed in [6]. 
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momenta (see Fig. 1) 

kt + k '+ 
p+ +p 

the skewness parameter 



p+ _|_ p'+ 

and Mandclstam t. Analogous definitions hold for the antibaryon vertex. The u ^ c 
transition GPDs are expected to exhibit a maximum which becomes more pronounced as 
the mass of the heavy quark increases. The position of the peak is approximately at 

xo = rric/M = 1 - e/M, (5) 

for not too large —t. The width of the peak shrinks with increasing mass. For standard 
values of the charm quark mass [12] 

m, = (1.27+J;;j)GeV, (6) 

Xq is about 0.5 - 0.6. The parameter e is the difference between the hadron and the heavy 
quark mass and has a value of about 0.6 — 1.0 GeV. In the formal limit of rric — > oo 
xq tends to 1 according to the heavy quark effective theory (HQET) [13]. This expected 
property of the transition GPDs parallels the theoretically expected and experimentally 
confirmed behavior of heavy quark fragmentation functions, in particular c A+ (see 
[14, 15]) and is also analogous to the behavior of light-cone wave functions (LCWFs) and 
distribution amplitudes for heavy baryons [16. 17]). The x dependence of the transition 
GPD is to be contrasted with wide-angle Compton scattering [2] and the zero-skewness 
GPD analysis [18] where the GPDs exhibit a behavior similar to the one expected for the 
u — > c transition only for —t^l GeV^. 

Some kinematical details of our process of interest are compiled in App. A. As for 
Compton scattering we work in a symmetric frame with the baryon momenta as specified 
in (88) and assume that the virtualities of all the partons the baryons are made of are 
smaller then A^ where A is a typical hadronic scale of order 1 GeV. Of course, the 
virtuahties of the charm and anticharm quark respect the relation | k''^ — ml |< A^. We 
further assume that all parton transverse momenta satisfy the restrictions: 

kl,/x,<A^ k'l,/x[<A\ (7) 

where, for the ease of legibility, we have not distinguished between the parton momenta 
from the upper and lower vertex. The bounds for the transverse momenta are assumed 
to hold in hadron frames in which the respective hadron moves along the 3-direction 
(the tilde characterizes momenta in the hadron-in frame, the hat those in the hadron-out 
frame). These hadron frames are reached by transverse boosts, cf. e.g. [19]. This is a 
transformation that leaves the plus component of any momentum vector a unchanged. It 
involves a parameter b'^ and a transverse vector and is defined by 



[a+,a ,a_L] 



, _ a_L • b_L a"'"b]_ a"*" 



(8) 
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The corresponding transformation for antibaryons requires transverse boosts that leave 
the minus components of the momenta unchanged. In the following we discuss the baryon 
vertex in detail. All the arguments and results given below do hold for the antibaryon 
vertex analogously. 

The active parton momenta are parameterized as 

h = [x^p+, KM±] , K = W^, K", K±] ■ (9) 
The individual momentum fractions xi and x[ are related to xi and the skewness by 

Transforming the parton momenta (9) to the respective hadron-in or hadron-out frame 
by the transverse boost (8), we find 

^1 = ...M± + xiAj2] , k[ = [x[p'+, ...,K±- x[Aj2] . (11) 

For convenience let us assume that there is only one spectator ^. The generalization to 
an arbitrary number of spectators is straightforward but tedious. As is characteristic of 
any parton approach the momenta of the active parton {ki and k'l) and the spectator {ks 
and kg) sum up to the parent baryon 's momentum. Hence, 

k±s = - kix , Xs = 1-xi, (12) 

and analogously for /c^. The spectator condition /c^ = A;^ leads to a;'^(l — ,^) = a;i(l+^) — 2,^ 
and k'^j^ = + kix. Transforming also the spectator momenta to the respective hadron 
frames and using the supposition (7), we can write 

2A^ > ~klJxs + k'lJA's = ^Ai/2 + -4-(ki± + A^/2)^ . (13) 

J- S i- Xi 

Each term in this inequality is positive. Therefore, each one is separately bounded by 2A^. 
The inequality for the second term on the right hand side of (7) compels kij_ ~ — Aj_/2 
for large rric- This result can be approximated by the more convenient form 

kij_ ~ — xiAj_/2 for rUc ^ oo (14) 

since the relevant values of xi are about xq, i.e. large. Analogously, one obtains k'j^j^ ~ 
x[A±/2. Approximating xi by the peak position xq, we find from the first term in (13) 

Al<4-^(l-a- (15) 

i- — Xq 

In the formal limit of rric —>■ oo where xq tends to 1, this inequality is respected for any 
fixed value of A^^. Even for a realistic value of about 0.6 for xq the range of allowed 



^This may be viewed as a quark-diquark configuration of the baryons. 
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values for is sufficiently large for our purpose of studying the process pp A+A~. 
Finally, the assumption of bounded virtualities, \ kj \< A'^ and | k'^ — ml |< A^, provides 
ki ~ xip" and an analogous result for k'l" . 

The above considerations lead to the following approximation of the momenta of the 
active partons (9) {k^^^ = 4'^" - 



ki 

K 

k2 



ki 

K 
k2 



1+ ^l^l 1 A ■ 



x|Ai 



5 ^2 ' 



2 



mg + a;fAi/4 _ 1 , 



(16) 



We introduce the subscripts 1 and 2 in order to distinguish the partons from the p — > 
A+ and the p — > A~ vertex. It is to be emphasized that the parton momenta in this 
approximation are on-shell and lie in the scattering plane formed by the baryon momenta. 



2.2 Calculation of the double handbag 

The amplitude for the hadron-hadron scattering process obviously reads, 

X (gV|r^:^„^(z2/2)*:,,,(-Z2/2)|gi^), (17) 

where we have used kf^ := {ki + k'j) /2, {ki — k[) = {p — p'), and {k2 — k'2) = {q — q')- It 
is understood that the proton emits a parton with momentum ki, helicity Ai and color 
ai (cf. Fig. 1). It undergoes a hard scattering with a corresponding antiparton (A;2, A2, 
02) emitted from the antiproton. The scattered partons, characterized by /c^, A'^, a'-^ and 
/c2, A2, 02' reabsorbed by the remainders of the proton or antiproton and form the 
A+ or A~, respectively. The labels af'' are spinor labels. For the ease of writing we will 
omit the color and spinor labels in the following wherever this can be done without loss of 
unambiguity. The active u and u quarks are treated as massless quarks. We have not to 
consider an antiquark contribution in the above baryon matrix elements as is necessary 
in Compton scattering [2]. The quarks in (2) are associated with the p — > A+ vertex, the 
antiquarks with the p A~ one. 

Since the parton-parton scattering is dominated by a large scale, the heavy quark 
mass, we can neglect the variation of transverse and minus (plus) components of ki and /c^ 
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in the hard scattering kernel H and replace the parton momenta by the on-shell vectors 
fc-'^ (16). With this replacement of the parton momenta the integrations over k^^^, k^^ 
and kj,!'^ ,i — 1,2 , can be performed exphcitly 

j d'kreikr) I ^^^'^^"/ d'k-e{kT-) I ^e^^t^-/(^i,.i,^„..) 

= J dkT^e{kT+) j ^e^'^'^^^ J dkre{kr) j ^e^'^'-^^fih,Zuk2,Z2). 

(18) 

We see that the relative distance of the fields in the matrix elements is forced on the 
light-cone, zi = [0, z^, Ox] and Z2 = [z^ , 0, 0_l]. After this the time-ordering of the fields 
can be dropped [20] . 

To show that the plus components indeed dominate the scattering we use the fact that 
the proton-parton amplitudes, described by the soft matrix elements, can be written as 
the amplitude for a (anti) proton with momentum p (g) emitting the active parton with 
momentum ki{k2) and a number of on-shell spectators times the corresponding conjugate 
amplitude for momenta p'{q') and k'i(k2)i summed over all spectator configurations; this 
corresponds to inserting a complete set of intermediate states between the quark and the 
antiquark fields in the matrix elements. 

Using the technique of transverse boosts, we can transform to a frame where ki has 
zero transverse and minus components: 

^1 = [A:i+,0,OJ. (19) 

Using the general result 

i(7-7+ + 7V) = l, (20) 
and the energy projector for quarks ^ 

u{h, XiMh, Ai) = ^1 • 7 = /c+7- , (21) 

Ai 

we can write 

*"(V2) = i(^-^+ + ^V)*"(W2) = ^5^[^^(^i,Ai)(iZ(^i,Ai)7+*"(W2)) 

^ Ai 

+ 7+«(fci,Ai)(«(^i,Ai)^"(;2i/2))] . (22) 

Due to the central assumption of restricted parton virtualities and transverse momenta 
(7) there are large plus components at the proton-parton vertex but no large kinematical 
invariant. Because of this the first term in (22) dominates over the second one and the 

^Wc use the Dirac spinors as defined in [10]. Tlie spinors refer to states with definite hght-cone 
heUcities [21]. 
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latter can be neglected. Since a transverse boost leaves the plus component of any vector 
unchanged, (22) holds in the overall symmetric frame too. Similar considerations for the 
outgoing partons run into complications with the quark mass. A transverse boost to a 
frame where k[ has no transverse component, leads to 



2 

mi 



(23) 



with a non-vanishing minus component due to the finite charm-quark mass. The energy 
projector reads 



2 

Y,<k[,KMK,K) = K-j + m, = ^;V + ^7+ + me (24) 
in that frame. For the charm field we can therefore write (with (7''')^ = 0) 

n-zi/2) = ^X][«(^;,a;)(«(^;,a;)7+*^(-;.i/2)) 

+ 7+{«(^;,A')(«(^l,A;)^^,(-;.i/2))-2me^'=(-V2)}] • (25) 

Using the plane-wave decomposition of the charm-quark field it can be shown that the 
term within curly brackets vanishes identically if it is applied to a Fock state in which 
the charm quark carries momentum k[. But these are just the states we need in the 
sequel to describe the absorption of the charm quark, that leaves the hard scattering 
process with momentum k[, by the A^. We therefore conclude again that the first term 
in (24) dominates over the one within square brackets so that the latter can be neglected. 
In the overall symmetric frame this term is small as well and we have 

= ^ E "(^1' ^1) (^(^'i' Kh^n-^im) ■ (26) 

Employing the light-cone spinors we can write for the helicity non-flip and flip config- 
uration at the p — > A+ vertex in the symmetric frame 



u{k[, —Xi)ia~^^u(ki, Ai) 



1 



2^/ktK 



^(2Al^)^■ , (27) 



where j = 1,2 and 



= ^^±y _ (28) 

Of course, analogous expressions hold for the antiparticle spinors (with an additional 
minus sign). The helicity projector for the light quarks reads 

u{ku XiHh: Ai) = h ■ 7(1 - 2Ai75)/2 . (29) 
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On the other hand, for the charm quark we have 

u{k[, X[)uik[, a;) = [^i • 7 + mc] (1 + 75-^1 • 7)/2 , (30) 
where the covariant spin vector Si reads 

Here, = [0, 1,0_l] is a unit vector. In the hmit rric — > the spin projector (30) turns 
into the usual hehcity projector. 

Now we can simphfy the products of quark fields occurring in (17). In fact, forming 
the product of the first terms in (22) and (25) and inserting (27) between the two terms 
involving 7+, one obtains after a little algebra 

K'S-^imKMm = ^^E(^'^(-V2)7+«(fci,A;)) 

X (u{h, Ai)7+M/"(^i/2))m„. (fc;, X[W,{ki, Ai) 



^y^Tf E {^'^(-^1/2) 7-^^^^^*"(^i/2)M(^i' ^i^aAku Ai) 

- i{-2Xiiy^%-zi/2)ia+^ ^^'f'^' ^^{zi/2)u^'fi[,^ . (32) 

The quark- mass and the transverse-momentum terms drop out since they always come 
with products of two 7"*" matrices. Analogously one finds for the p ^ A~ vertex 

X |^"(V2) 7"^^|^^1-^2/2)i;a,(^2, A2)t;a^(^^, A2) 

- z(-2A2z)^'^"(;^2/2) la-^ ^ + ^^^^' ^^(-;^2/2)t;.,(fe2, X2W'M, -Aa)) .(33) 



Combining now (17), (18), (32) and (33) and using A2 = — Ai, a property of the subprocess 
amplitude which follows from the fact that massless quarks do not flip their helicities in 
interactions with photons and/or gluons, one arrives at 



9 



27r A^ktk'+k^k'^- 



Ai 

X (gV I ^r«(^2/2)7-^-^^^'=(-^2/2) I gz.)i/A,-A„A.-A,(^-i,fc2) 

- (pV I ^^(-^i/2)z<7+^i±|^*"(^i/2) I p^i) 

X (gV I ^"(^2/2)za-^' ^~^^^^^ ^--(-Z2/2) I qiy) H^x,M,M-M{h: 

- l{2X,iy {p' 1^' I ^rC(_^^/2)^+l±|^^«(^^/2) I p;,) 

X (gV I ^-(52/2)a-^'^^^|^*^(-Z2/2) | qv) Hx,M,M-M{h, h) 

- z(-2Aiz)^(pV I r(-5,/2)<j+^i±|^*-(zi/2) I P//) 

X {qV I fr«(^,/2)^-i±|^^^(-^2/2) I gi.)i/_A,-A„Ai-A,(^i,^2)} . (34) 

The subprocess amplitudes are defined by 

Hx',x',MxAki, k2) = v{k2, \2)u{k[, X[)H{U;\ k^2^)v{k'„ X'Mh. Ai) , (35) 

where the hard scattering kernel H(k^^\ ^2'*) has to be evaluated for parton momenta as 
given in (16). 

2.3 Peaking approximation, GPDs and form factors 

According to our discussion in Sect. 2.1 we set k^^^^ = x^pp^'^^, k^' = x'^q'^'^'. In 
order to produce a cc pair the subprocess Mandelstam variable, s — X1X2S, must be 
larger than 4m^. For the sake of the argument consider a symmetric configuration of 
p —>■ A+ and p A~ which, with regard to the expected behavior of the GPDs, namely 
a marked narrow peak centered around Xq, is an appropriate assumption. We then have 
the kinematical requirement 

M , , 

X > 2xo-= . (36) 

As the comparison with (94) in the appendix reveals, the inequality (36) requires x to be 
larger than the skewness in the forward hemisphere. With the help of (10) one sees that 
the inequality (36) also implies x > 2(1 + C,)xqM / ^ — and, hence, x > ^. Consequently, 
the ERBL region a; < does not contribute to the process pp — > A+A^. For large s the 
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lower bound of x becomes x > 2xoM/y/s, i.e. it tends to zero for s — > oo as does ^, see 
Eq. (94). 

Given the expected shape of the GPDs it is clear that only regions close to the peak 
position contribute to any degree of significance. For the hard partonic subprocess we 
therefore employ the peaking approximation, i.e. for the momenta (16) we use Xi ~ a;2 — 
Xq. From the parton momenta (16) it then follows that the Mandelstam variables for the 
subprocess are proportional to those of the full process 

s ~ XqS , u ~ xlu , £ ~ t . (37) 

The sum of the three subprocess Mandelstam variables is 

s + u + tc^ 2ml ■ (38) 

Thus, momentum conservation for the parton kinematics holds. The use of the peaking 
approximation is quite appealing since it nicely matches the subprocess kinoiiiatics to the 
one for the hadronic process up to corrections of order m^/s and ,^(1 — xq). The accuracy 
of the peaking approximation increases with increasing s at fixed mass of the heavy quark 
and with increasing rric for fixed values of m^/ s. For the physical charm system the peaking 
approximation holds, for instance, for (38) or the cosine of the c.m.s. scattering angle 6 
with an accuracy better than 10(20)% for s = 30(25) GeV^ in the forward hemisphere. 
This means that cos 9 cither evaluated from s and t or from s and t only differs by less 
than the quoted accuracy. The use of the peaking approximation therefore implies that 
our approach can only be applied to values of s well above the kinematical threshold. The 
natural requirement that a substantial part of the peak region of the GPDs contributes 
to the form factors necessitates $, <^ xq which in turn implies not too small values of s. 
Numerical checks reveal that this requirement is satisfied provided s^23 GcV^. 

The Fourier transforms of the soft hadronic matrix elements ^ occurring in (34) define 
the GPDs as we discuss in App. B (104), (119) and (129). Using 



/ 



.av+^(^r) _ f dx. 



dkf'^ ' ' = / , ' , (39) 



and an analogous expression for the integral over as well as the covariant decompo- 
sitions (108), (122) and (129), one sees that the integrations lead to 1/xi moments of the 
GPDs. For instance, for the vector matrix element we arrive at 



- Rv{^^ tMp'^ ^I'hMp, + M^, t)u{p', ^i')^-r-^u{p, /.) . (40) 



^In order to have manifestly charge conjugation symmetric GPDs (see App. B) we have to add the 
corresponding antiparticle contribution to the bilocal matrix elements of the field operators in (34) (see 
definitions (103) and (104) and analogous ones for the other GPDs) which can be done since these 
contributions are zero because of the neglect of intrinsic charm. 
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The form factors in (40) are defined by 




dxi 



=Fi{xi,i,t) , 



(41) 



where and Ri are generic for any of the u ^ cor u c transition GPDs H, E, H, E, Ht, 
■ ■ ■, and the associated form factors Ry, Rt, Ra, Rp, St, ■ ' We apply here the notation 
proposed in [22]. Due to the kinematical requirement x > ^ the lower limit of the in- 
tegration is taken to be ^ instead of zero. This leads to a mild ^ (or s) dependence of 
the form factors. For s oo the ^ dependence of the form factors disappears since in 
this limit ^ tends to zero like (M^ — m?)/s, as can be seen from (94). The form factors 
Ri stand for linear combinations of more commonly used form factors, as can be seen 
from Eqs. (Ill) and (125) in the Appendix. Only such linear combinations occur in the 
process amplitudes. Due to the time reversal invariance the transition form factors are 
real valued. 

It may seem plausible to choose the peaking approximation for the factor 1 / \Jxl — ^ 
in Eqs. (34) and (39) too and to use for the form factors the expression 



For GPDs concentrated in a narrow range of large centered around Xq, the difference 
between (41) and (42) is expected to be small. We will return to this issue when we 
present our model GPDs. 

Working out the various hadronic covariants occurring in (40) and the other baryon 
matrix elements in the symmetric frame, we find the following results for quark helicity 
non-flip (conventions - helicities of the outgoing (anti)baryons appear always flrst; x is 
either x\ or ^2) 




(42) 






MRp , 




(43) 



M + m 

(— l)'^"^/^, S = (— )1 for (anti) baryons. 



^ (l + OM-(l-Om 



RveS — Rv — Rt 



1-e 



■M. 



(44) 



Here, e — 
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For quark helicity flip we define the combinations (A^ refers to tlie lielicity of tlie 
emitted parton) 



n^Tcu 
njTm 



n_/Tcu o \ -ijTcu 
'TjTcu r) \ a_/Tca 



(45) 



of helicity-fiip GPDs [23] which are defined in (129) and (130). These combinations make 
the hehcity projection exphcit since they correspond to the Dirac structure i/2a^^{l ± 
2A175) (see (132)). Their moments read 



1: 



dx 



Tcu 



/-T/J CU \ 

2 Ai/i/i V ^Ai/i/i/ 



5- 



P 



+ 



I (M + m) 



dx 



njTcu 
2 '^M-HH 



+ [1 + 2X,Se^] Sv, - [e + 2Aife] Sy, } , 

^^{il-e){Se + 2X,)Sr + Se^Ss 



^5 



(5s + 2Xi)MkSv,-Sv, 



(46) 



and represent the quark helicity- flip matrix elements occurring in (34). The form factors 
are defined in (41). 

The matrix elements (43) and (46) vanish for t ^ to a,t least as 



OC \/to —t 



|Ai-//'-Ai+/u| 



(47) 



This behavior reflects the fact that the mismatch of the quark and hadron helicities in 
these matrix elements is to be compensated by a corresponding number of units of orbital 
angular momenta in order to ensure angular momentum conservation. Closer inspection 
of the matrix elements for the various helicity configurations reveals that only Ry, Ra and 
St (corresponding to the GPDs H, H and Ht) do not involve orbital angular momentum 
and therefore do not vanish for t ^ to- The form factors Rt, Rp, Syi and require 
one unit of orbital angular momentum while Ss (corresponding to the GPD Ht) requires 
even two. 



3 Amplitudes 

3.1 Subprocess amplitudes 

The amplitudes Hy y ,AiA2 for the subprocess (2) calculated to lowest order of perturbative 
QCD, read 

^ 8TTas{s) ui — it 

~ ~s ^1 - Aml/s ' 
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s 



= (48) 

Since we consider the u quark as a massless particle all amplitudes with the same helicities 
of the u and u quark are zero. Helicity amplitudes other than those given in (48) are 
obtained by applying the relation 

-AiAi = -(-l)^'"^^-f^A'iA^,Ai-Ai , (49) 

which follows from parity invariance. With the help of the peaking approximation we can 
express the subprocess amplitudes in terms of the kinematical variables of the full process 

= — 87rQ;s(a;o>s) cos^ , 
= 87ra^(xos) sin^^/2, 

= 87ras{xQS) —^sinO . (50) 

The quahty of the peaking approximation has been discussed in Sect. 2.3. The QCD 
coupling constant as is evaluated for four flavors with Aqcd = 0.24 GeV. 
The cross section for the subprocess reads 

da 1 Ajn, . ._,s 

where A^^ is defined in (92) and 



1 



^0 = I I Hxix',M-x, I' = (47ra.)' [l + cos^ 6 + ^ sin^ . (52) 

To lowest order of QCD the subprocess amplitudes are real. As a consequence any 
transversal polarization of the involved quarks is zero since it is given by the imaginary 
part of a flip-non-flip interference. To higher order of QCD this will change. The sub- 
process amplitudes become complex [24] and a c-quark polarization of order a^mc/ 
is obtained. The polarization of the u quark is zero to all orders since it is treated as 
massless. 



AM 



2 



3.2 The hadronic amplitudes 

Inserting the definitions (104), (119), (129) as well as (130) into (34), we obtain for the 
helicity amplitudes of the process pp — > A+A~ 



Cf\ 1 f dxi dx2 



Ncj 8(p+)2 J y^F^ 
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Here, the ratio Cp/N^ is a color factor (=4/9). In principle, one may express the 
amplitudes in terms of the form factors using (43) and (46). However, this will lead to 
a lengthy expression in terms of the eight form factors. We refrain from doing this. The 
following hierarchy of these form factors is to be expected 



St Ri 



Rt 



s- 



VI 



V2 



(54) 



This expectation is based on the fact that the supposedly suppressed form factors involve 
orbital angular momentum of the quarks making up the baryons. It is also supported by 
the model calculations for the form factors which we will present in Sect. 4. The analogous 
hierarchy has also been found for the form factors parameterizing the elastic p — > p matrix 
elements which, for instance, occur in wide-angle Compton scattering [2]. The reason for 
it is that according to theoretical expectations and phenomenological experience, the 
valence Fock states of many hadrons are dominated by parton configurations with zero 
orbital angular momentum. Given the hierarchy (54) and the smallness of skewness (94), 
we may neglect all terms oc ^ in (43) and (46). This way the form factor Rp decouples 
completely. For quark helicity flip wc have additional suppression from the subprocess 
amplitudes (48) in the forward hemisphere. As can be seen from (53) configurations with 
quark helicity flip at both hadronic matrix elements come with the amplitude +_ and 
are therefore suppressed by an additional factor — t. Interference terms between flip and 
non-flip matrix elements are proportional to -tip and are suppressed by a/£o — irric/s. 
Due to these properties one may neglect all form factors related to the helicity-flip GPDs 
with the exception of St- The model GPDs and form factors which we will present in 
Sect. 4 support this assumption. 

For later use we provide the explicit expressions for the helicity amplitudes by taking 
into account only the most important form factors Ry, Ra, St- Inserting (43) and (46) 
into (53) we obtain the following set of helicity amplitudes for pp A+Aj 

Cf 
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+,+- 
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Ry^ Ra. 



C-2 
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Ri 



Ra 



(55) 
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The remaining amplitudes are fixed by parity invariance 

M_^,_,,,_^_. = M^v',M. ■ (56) 

One can easily convince oneself that the combination Ry + Ra parameterizes the soft 
hadronic matrix element for the situation where the proton with a given helicity emits a 
u quark of the same helicity and subsequently a c quark with the same helicity as the u 
quark is re-absorbcd forming a with the same helicity as the proton. On the other 
hand, the combination Rv — Ra stands for the situation where the proton and the A+ have 
still the same helicity but the quarks have now a helicity opposite to that of their parent 
baryons at one vertex either at the p — > Ac or the p — > Ac one. The other amplitudes take 
into account configurations with either one or two p — > A+ helicity flips. 

Strictly speaking, our helicities are light-cone helicities. A standard c.m.s. helicity 
amplitude is given by the light-cone helicity amplitude with the same helicities as the 
c.m.s. one plus an admixture of other light-cone helicity amplitudes 

(57) 

The strength of the admixture is controlled by the parameter (3 which reads for our 
kinematics [23] 

^ _ 2M ^tu - tpuo ^^g^ 

(1 + Am)^ AmV-su - u - ^/touo ' 

where the proton mass is neglected. Near the forward direction (3 is small and its effect can 
be ignored, but for larger scattering angles it is to be taken into account in the calculation 
of spin effects. For <^ s (58) simplifies to 

P^^^y^. (59) 

In this approximation /3 has been used in wide-angle Compton scattering [25]. 



4 Modelling the GPDs 

We will model the GPDs by overlaps of LCWFs for the A+ and the proton. Only the 
valence Fock states of the baryons are taken into account by us. For the A+ this is 
probably a reliable approximation. For the proton, on the other hand, higher Fock state 
contributions are likely important but the required overlap with the Ac wave function 
projects out only its valence Fock state. 

The overlap formahsm for GPDs in terms of LCWFs has been developed in [26, 27]. 
For the process of interest there is a little comphcation with it. Due to the unequal mass 
kinematics p'^ becomes negative for large scattering angles in the symmetric frame. The 
spectator partons, in this case, move in the direction opposite to their parent hadron. In 
order to apply the overlap mechanism in this situation one would have to boost first into 
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a frame in which pg is positive. Although this can be done in principle, it is fortunately 
unnecessary since the pp — > AcAc cross section is strongly peaked in the forward direction 
for not too small values of s, as will turn out below. Thus, it suffices to work out the form 
factors and the cross section only in the forward hemisphere where p'^ is positive. 

In the ERBL region where a; < ^, meson poles, here D and D*, may contribute to 
the GPDs. However, as we discussed in Sect. 2.3, this region does not contribute to 
our process since small x are insufficient for the generation of the cc pair. One may 
also think of Regge exchanges in general, in addition to the handbag contribution. The 
most prominent Regge pole that can contribute to our process is again the D*(2007). 
Taking a standard slope for the corresponding Regge trajectory one finds an intercept of 
a{Q) = —2.63 for it. Since the cross section, say at t = 0, behaves as da/dt s2(a(o)-i) 
for a Regge exchange its contribution is strongly suppressed as compared to the handbag 
contribution. Thus, it seems save to neglect pole terms. 

In order to calculate the GPDs within the overlap formalism we need the LCWFs 
of the proton and the A,.. For the proton a number of such wave functions can be 
found in the literature, e.g. [28, 29, 30]. The proton wave function advocated for by 
Bolz and KroU (BK) [28] has been determined by fits to the Dirac form factors of the 
proton and the neutron at large momentum transfer, to the valence quark distributions 
at large x and to the decay width of the process J/^ — >■ pp. It has already been applied 
in an overlap calculation of the proton GPDs at large —t which form the soft-physics 
input for an analysis of wide-angle Compton scattering within the handbag approach 
[2, 25]. The results of this calculation are found to be in fair agreement with the data 
on Compton scattering [31]. Recently the BK wave function [28] has been supported by 
results from lattice QCD. The first moments of the proton distribution amplitude, i.e. 
the LCWF integrated over the quark transverse momenta, that have been calculated in 
[32] , are in good agreement with those evaluated from the BK wave function. It should be 
noted, however, that higher moments tend to weaken the asymmetry in the distribution 
amplitude that one would expect from the first moments only [33] . The nucleon DA from 
the lattice could thus be closer to the asymptotic form of a spin 1/2 baryon DA than 
the nucleon DA resulting from the BK wave function. We think, however, that this is of 
minor importance for our pTirpose of modelling the m — > c transition GPDs and we will 
use the BK LCWF in the following. 

The valence Fock state of the proton can be expressed as 



p, +; uud) = / [dxUd'^kj]^ {^123^+-+ + *2i3A^-++ - (^132 + *23i)A^++-} , (60) 



with the plane wave exponentials and the color wave function omitted. The integral 
measures are defined by 
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(61) 
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KK (66), (67) 


2117 


0.75 


0.90 


460 


p = 2.0 


BB (66), (68) 


3477 


0.75 


0.90 


431 


p = 2.0 


(77) 


169.5 


1.1 


0.90 


317 


diquark 



Table 1: The parameters of the hght-cone wave functions used for the calculation of the 
u ^ c transition GPDs. For the Ac the quoted r.m.s. transverse momentum is that of the 
heavy quark. 



A three-quark state is given by 

A^AiAaAa = , | Xi, k_Li, Ai) | X2, k_L2, A2) I ci;x3,k_L3, A3) . (62) 

1^X1X2X3 

The single particle states satisfy standard normalization. Eq. (60) is the most general 
ansatz for the L3 = projection of the three-quark proton wave function. Prom the 
permutation symmetry between the two w-quarks and from the requirement that the 
three quarks have to be coupled in an isospin 1/2 state it follows that there is only one 
independent scalar wave function [34] which according to [28] is parameterized as 



'^l^{Xi, k_Lj) = *(xi, X2, X2,] k_Li, k_L2, k_L3) = Np{l + ZXi) CXp I ^ 

Xi 



(63) 



The Gaussian factor is typical for such light-cone wave functions. The distribution am- 
plitude that corresponds to the wave function (63) is 

^fi = 60x1X20:3(1 + 3x1) . (64) 

As usual it is normalized in such a way that J[dx]2^ — 1. The distribution amplitude 
(64) is displayed in Fig. 2. It is rather similar to a distribution amplitude obtained from 
QCD sum rules [35]. The parameters Np and of the wave function (63) have been 
determined in [28]. They are quoted in Tab. 1 together with the r.m.s. value of the quark 
transverse momentum and the probability of the proton's valence Fock state evaluated 
from this wave function. 

According to the HQET [13] the major part of the Ac state comes from the config- 
uration where the c quark carries the helicity of its parent baryon while the two light 
quarks are coupled in a spin-isospin-zero state. In the formal limit of ttIc 00 there 
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Xi 



Figure 2: The distribution amplitudes (64) (dotted), (67) (solid) and the one correspond- 
ing to the wave function (68) (dashed line) versus Xi at X2 = x-^. (colors online) 

should be no admixture from configurations where the helicity of the c quark is opposite 
to that of the A^. However, for a c-quark mass of 1.27 GeV which is large enough to 
justify the application of HQET ideas but too small in order to suppress corrections to 
the HQET decisively, a small admixture from a c-quark with helicity opposite to that of 
the Ac cannot be excluded. We therefore write the valence Fock state of the Ac in the 
form 

|A+,+;cMrf) = y"[rfx]3[rfV]3{(A^;+--A<;-+)+p(a;2-X3)A<%+}^A(a;i,kxi), (65) 

where are the three-quark states defined in (62) with the first u quark being 

replaced by the c quark. The scalar wave function, \I/a, must be symmetric under the 
interchange of the momentum fractions of the two light quarks. The factor X2 — in (65) 
takes care of the isospin of the p term [36] . It also guarantees the vanishing of the p-term 
in the formal limit iric — >^ oo as is required by the HQET. 

We know from the LEP experiments [37, 38, 39] that the polarization of the h quark 
generated in the process e+e^ Z ^ bb (P^ = —0.94) is not fully transferred to the 
Af, baryon as is predicted by the HQET. Rather the A;, polarization only amounts to 
—0.43 ± 0.25 as an average of the three experimental results. A substantial fraction of 
the observed Ab baryons likely comes from the decay of the Sf, and which leads to 
a strong depolarization of the At,. Whether this explains the observed A;, polarization 
fully or only partially is unknown. Hence, the possibility that there is an admixture from 
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configurations where the and the b quark have opposite helicities cannot be ruled out 
at present. It is plausible to expect that for the lighter c-quark system this configuration 
is even more important. We therefore allow for it and take a value of 2.0 for the mixing 
parameter p. The configuration with opposite helicity then contributes about 10% to the 
probability of the state (65). 

We parameterize the scalar wave function in (65) as 



A^A exp [ - f{xi)] exp - ai ^ 



2-, 



(66) 



The function f{xi) has to generate the expected pronounced peak at Xi ~ Xq as defined in 
(5). Current model wave functions for the Ac possess this property [5, 16, 17, 40, 41, 42]. 
For f{xi) we adopt a slightly modified version of the one given in [16], 

fMx^) ^alM'^f^. (67) 



Xi{l - Xi) 

The LCWF (66) is a suitable adaptation of the one proposed in [41] which has been 
obtained by transforming a harmonic oscillator wave function to the light cone [42]. 

As an alternative we will also make use of a result from QCD sum rules obtained for 
the Ab baryon [17] which we adapt to our case of a charm baryon, namely 



fBB{xi) = aAM(l - Xi) . 



(68) 



Up to corrections of order (M — mc)/M the two functions have the same Xi dependence 
for large rric. The parameters of the two choices for the Ac wave function are quoted in 
Tab. 1. We assume Oj, = oa and a probabihty of 0.9. The normalization A'a is fixed 
from these parameters. A large value of the Ac valence Fock state probabihty is to be 

expected [41]. 

The distribution amplitude corresponding to (66) reads 



$A oc X1X2X3 exp [ - f{xi)] . 



(69) 



The two Ac distributions amplitudes are compared to each other and to the proton distri- 
bution amphtude in Fig. 2. The variant (68) is broader than (67) and is very similar to 
the proton distribution amplitude. Since the c quark is relatively light the peak positions 
of the Ac and the proton distribution amplitudes do not differ much. For the case of the 
Ab, for instance, the peak positions would differ markedly. 

The GPDs are obtained from the following overlap integrals 



H^^{x,i,t) 



1 



X 



^9 I [c?a;]3[d^k_L]35(x - Xi)*A(^i,^±j) 

^123(5*, ku) + ^132(^1, k±_i) + ^^231(^1, ku) 



AH^{x,tt) = 

X 



p 



[dx]3[d'^k±]3S{x - xi)*A(^i, k±i) 



1_^2 

(£2 - £3)^213(^1, k±i) , 



(70) 
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and the GPDs H and H are given by the combinations 

All other GPDs are zero. They require parton orbital angular momentum, a property 
that our simple wave functions for the Ac do not possess. At least for large rric its struc- 
ture follows from the HQET. The proton wave functions might have a more complicated 
structure involving, for instance, components with non-zero orbital angular momentum. 
This would still lead to the same result for the GPDs since the overlap with the Ac LCWF 
does project out only L = components of the proton wave function. 

The arguments of the light-cone wave functions are related to the average parton 
momenta, over which it is to be integrated, by [26] 



xi + i 7 J 1 — xi A 



Xj — — ^ , k±j = k±j H ^ — — , for j = 2, 3, 



Xi = 



k±i 


= k±i - 


k±j 


= k±j + 


k±i 


= ^±1 + 




= k±j - 



1 + ' ' 1 + ^2' 

xi-^ r r 1-XiAj 



1 - e ' 1-^2 



% - k±j ^ ~k^j - -^^^ , for J = 2, 3. (72) 

Working out these integrals for the wave functions (63) and (66) we obtain the GPD 
3 NaNp 1 x^-C^ 



HI 



cu 



4 (47r)M-ea^(l + + 4(1-0 
{l-xf 1 + 2^ + S 



X 



X exp 



1 + e a2(l + 0^(x -0 + 4(1 -0^(^ + 

(1-x) aWpAl 

a2(l + 0^(^-0 + 4(1 -0^(^ + 0- 



f(^-^\\ exp r -(l-xjaja^Aj n 

^ U _ [alii + mx - £) + alii - mx + £) J ' 



and for the p-dependent part 

15 1 - C 1 + 2C + S ^ ^ ' 

The GPD Ht is shown in Fig. 3. One observes that the GPD exhibits the anticipated 
pronounced peak at a; ~ Xq. With increasing (or -t', see Eq. (98)) the position of the 
peak is slightly shifted towards larger values of x. This property of the m — > c transition 
GPDs is similar to that of the proton valence-quark GPDs which have been determined 
from an analysis of the nucleon form factors through the familiar GPD sum rules in Ref. 
[18]. The only difference is that the it — > c transition GPDs possess this property already 
at t' = as a consequence of the intrinsic scale Am^^. The proton GPDs exhibit a peak 
only at large values of —t and, therefore, the handbag approach can be applied to (real) 
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Figure 3: The u c transition GPD Ht obtained from the overlap of the LCWFs (63) 
and (66) (for the variants (67) (left) and (68) (right)) versus x at s = 30 GeV^ and A^^ = 
0,2.0,4.0 GeV^ corresponding to t' = 0, -2.17, -4.41 GeV^ (^ = 0.112,0.123,0.139) 
shown as solid, dashed and dotted lines, respectively (colors online). 

Compton scattering only at wide angles. The two variants of the mass exponential entail 
striking differences in shape of the resulting GPDs displayed in Fig. 3. The GPD evaluated 
from (68) is much broader than the other one. The maxima of the GPD evaluated from 
(68) are sited at lower values of x for small A^^ but at slightly larger ones for large A\ 
than is obtained from (67). The differences between the two GPDs have a bearing on the 
form factors and predicted cross sections as we are going to discuss in the following. 

With the GPDs at hand the transition form factors can be evaluated according to 
(41). They exhibit a structure analogously to (71) 

Rv{^, t) = St{^, t) - AR{^, t) , i?^(e, t) = ST{i, t) + Ai?(e, t) . (75) 

In Fig. 4 the form factors 5*^ and A/2, scaled by t' for the ease of graphical representation, 
are shown at s = 30 GeV^. The differences between the models (67) and (68) are partic- 
ularly substantial at small A^ or —t'. Also displayed in Fig. 4 is the approximation (42). 
As can be observed it is very close to the full form factor. This approximate coincidence 
nicely demonstrates the internal consistency of the peaking approximation. As the figure 
also reveals, Ai? is tiny despite the fact that the p term in (65) contributes as much as 
10% to the probability of the A+ state. The strong suppression of Ai? is a consequence 
of the tiny overlap of the p term in the Ac state with the proton wave function which 
manifests itself in the factor (1 — xY in (74). Even by a doubling of p, which leads to 
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Figure 4: The transition form factor St scaled by t' versus t' at s = 30 GeV^ evaluated 
from (73) using (67). The band indicates the parametric uncertainties of the form factor. 
The dashed line represents the the result obtained with (68) instead of (67) and the dotted 
line the approximation (42). The lower dashed line is 10 times the corresponding scaled 
form factor AR (colors online). 

the implausibly large contribution of 30% to the probability of the Ac valence Fock state, 
we obtain | AR |^| St \- Thus, AR is so small that it has no direct bearing on the 
observables for realistic values of p. 

Let us now turn to the issue of the error assessment. It is true that the independent 
parameters of the LCWF are the normalization and the transverse size parameter, 
but we do not have any idea about their uncertainties in contrast to the probability and 
the r.m.s. transverse momentum. We, therefore, assume errors for the latter quantities 
and vary A^a and a\ as well as rric (within the PDG limits (6)) in such a way that the 
errors of Pa and (A;^)^/^ are covered. For the latter quantity we assume an error of ±10% 
which is in fair agreement with an estimate of it for heavy baryons made in [40] . For the 
probability of the Ac valence Fock state we allow for the maximum value of 1 and take as 
its lowest value 0.7. The described estimate of the parametric error is shown for the form 
factor St in Fig. 4. We do not take into account the uncertainties of the proton wave 
function. They are small compared to that of the Ac wave function since the proton wave 
function is determined from detailed fits to nucleon form factors and parton distributions 
[28]. The uncertainty of the form factor evaluated from the mass exponential (68) is of 
about the same size as for the variant (67). 
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It is popular to model the baryons as a quark-diquark state. For the case of the Ac this 
appears rather natural with regard to the HQET and indeed many examples of quark- 
diquark wave functions for the Ac can be found in the literature [5, 16, 40, 43, 44, 45]. 
In all these cases only a spin-isospin scalar diquark, S[ud\, is used, an assumption that is 
in accordance with the large-mc limit of the HQET. In view of our experience with the 
three-quark models we omit a possible admixture of other diquarks for simplicity. The 
proton may have a more complicated spectator structure with scalar and axial-vector 
diquarks and perhaps excitations of higher states involving angular momentum [46] . But 
as we mentioned in the discussion of the 3-quark models, the overlap with the A^ LCWF 
only projects out the corresponding component with the scalar diquark from the proton 
LCWF. In order to render possible a comparison with the results presented in [5] we 
also employ the quark-diquark model here. For the ease of comparison we use the same 
parameters as in [5] although the LCWFs are shghtly modified in order to take care of 
recent developments. 

According to [5] we write the wave function of a proton in an u±S[ud] Fock state as 



*p(qD) = Np^qD){l " x) CXp [ 



x{l - x)- 



(76) 



The parameters of this wave function are quoted in Tab. 1. The probabihty of 0.5 for 

a proton may appear rather large. However, one may argue that a diquark is a bound 
state of two quarks and embodies therefore also gluons and sea quarks. Hence, the quark- 
diquark wave function effectively takes into account higher Fock states and has therefore a 
larger probability than one would expect for a 3-quark valence-Fock-state wave function. 
In view of this fact a larger transverse size of the quark-diquark state appears plausible, 
too. 

For the LCWF of a A+ in a c±S[udi state we use 
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x) exp 



2 y>^2l^-^0j 
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x{l — x) 



(77) 



x(l — x) - 

The parameters are again listed in Tab. 1. 

The overlap of the two wave functions (76) and (77) provides the following result for 
H — H — Hj' 
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(78) 



The other GPDs are zero. The quark-diquark GPD has a similar shape as the 3-quark 
GPDs. Only its size is larger by about a factor of 1.6. The form factors Ry = Ra = St 
in the quark-diquark model also behave similar to those in the 3-quark model. 
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5 Observables 



5.1 Cross sections 

The differential cross section for tiie process pp A^A~ reads 

da 1 . . da 11 Am .-^^ 

^ = ^.A.AM^ = ^-^ao, (79) 

where 

^0 = ■ (80) 

The cross section can readily be calculated from the amplitudes (55) using the subprocess 
amplitudes given in (48) and appropriate form factors. Particularly simple are those GPD 
models in which the Ac only consists of the configuration where the two light quarks are 
coupled in a spin-isospin-zero state (bound in a di quark or not) as is predicted by the 
HQET. They do not involve quark orbital angular momentum and the parameter p is 
zero. As we discussed in the preceding section a non-zero value of p is possible but cannot 
be large as estimated from measurements of the A;, polarization at LEP [37, 38, 39]. In 
the cases we examine, the p term does only lead to tiny differences in the form factors. 
Practically, there is only one independent form factor 

St{U) Rv{^,t) Ra{U) . (81) 

All other form factors are zero. In a situation where (81) holds, the differential cross 
section strongly simplifies and becomes proportional to the subprocess one 

%i,,^KK) - (g\:-,TO^Mi^ ,s.) 

Results for the differential cross sections are shown in Fig. 5. The differential cross 
section is sharply forward peaked. This effect becomes more pronounced for higher en- 
ergies. It arises from the dominance of the amplitude A4+-,+- (oc cos^^/2) near the 
forward direction and is natural given that we consider an annihilation reaction at high 
energies. The effect is pronounced by the behavior of the form factors which obey (81) and 
decrease monotonically with increasing —t'. The cross section increases with s for small 
scattering angles near the forward direction. This is a transient effect which disappears at 
large s and da/dfl ultimately show the usual 1/s decrease since the form factors become 
independent on ^ in that region as we mentioned in Sect. 2.3. The various models we 
employ behave similar in that respect. The quark-diquark model predicts a particularly 
steeply falling cross section. This feature is only to be expected given the small values of 
the r.m.s. transverse momenta in that model or the associated large transverse sizes (see 
Tab. 1). It can also be seen in Fig. 5 that the use of the mass exponential (68) instead 
of (67) provides markedly larger differential cross sections for small-angle scattering. 

The integrated cross section is displayed in Fig. 6. Its magnitude is of order nb. 
As the behavior of the differential cross section indicates (see Fig. 5) the largest cross 
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Figure 5: Predictions for the differential cross section versus cos 9. Left: s = 30 GeV^. 
Solid line represents the result evaluated from (73) with the mass exponential (67) (with 
error band). The dashed line is obtained with (68) and the dotted one is the results from 
the quark-diquark model. Right: The differential cross section at s = 23 , 30 and 50 GeV^ 
shown as dashed, solid and dotted line, respectively. The cross sections are evaluated 
from (73) with the mass exponential (67) (colors online). 

section is obtained with the mass exponential (68). For comparison we note that in [5] the 
pp A+A~ cross section was estimated to amount to about 3(10) nb at s = 30(50) GeV^. 
This is in reasonable agreement with the predictions from the quark-diquark model as well 
as the other ones presented here. 

In Sect. 4 we have already discussed the model uncertainties of the form factors. These 
uncertainties enter directly into the errors of the cross section (note that it is related to 
the fourth power of the form factor). The errors are shown as bands in Figs. 5 and 6. The 
width of the bands correspond to a probability of the relevant p-independent term in the 
Ac state (65) varying between 0.9 and 0.6 and to an error of ±10% for the r.m.s. charm- 
quark transverse momentum. All the models we exploit have similar parametric errors. 
Larger than these parametric errors is the difference between the predictions obtained 
from the different mass exponentials which is rather representative for the uncertainties 
of our predictions. 

5.2 Spin dependence 

For spin-dependent observables we have to take care of the fact that experiments provide 
typically c.m.s. observables which are expressed in terms of c.m.s. helicity amplitudes ($). 
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Figure 6: The integrated cross section (cos 6' > 0) versus s. For notation we refer to the 
left hand side of Fig. 5. 



We therefore have to transform our hght-cone amphtudes (53), (55) to the c.m.s. ones 
(57) with the help of the parameter (3 given in (58). We note that for observables that 
only refer to the helicities of the initial state baryons this transform has no effect up to 
corrections of order mj^. We also remark that the unpolarized cross section can be 
worked out in either basis, it is independent on the transform (57). 

While all single-spin asymmetries are zero to lowest order of perturbative QCD, there 
are many non-zero spin correlation parameter in our approach. We employ the conventions 
for the notation of such observables advocated for in [47]. The directions of the spin of the 
particles are denoted by L . . . helicity, iV . . . normal to scattering plane, S . . . 'sideways' 
spin direction within scattering plane but orthogonal to direction of momentum. 



Li = VilVi 



N 



N X L,- 



The corresponding spin states for and S type polarizations read 



IT(i)) 



1 

7! 



+) + (-)^ I -) 



1 

71 



+) + (-)!-) 



(83) 



(84) 



respectively. The spin-spin correlation parameters are denoted by (z, j = L, A^, S\. 

- A 



— C 



initial state spin correlations 
final state spin correlations 

polarization transfer from p(p) to A+(A^) — Dij(Dij), 

polarization correlation between p(p) and A~(A+) — Kij{Kij). 

A spin correlation parameter is defined by 

a{ij) + a{-i - j) - cr{i - j) - a{-ij] 



O 



1-3 



a{ij) + a{-i - j) + a{i - j) + a{-ij) 



(85) 



27 



where a{ij) stands for the differential cross section for the scattering with two particles 
in given spin states. For all GPD models that approximately possess the property (81), 
the spin correlations simply read ^ 

0{pp A+A+) = 0{uu^ cc) . (86) 

The form factors and, hence, their uncertainties, cancel out in the spin correlations. Thus, 
all our GPD models lead to the same set of predictions for these observables. In this sense 
our predictions for spin correlations are model independent. 

Below we list a set of examples of spin correlations obtained in the handbag approach 
from form factors which satisfy the relation (81). If the parameter /3 is only taken into 
account up to linear order, which is justified for not too large scattering angles, the 
corresponding analytical expressions become rather simple: 

All = -1, 

^ _ 1 + cos^ 9 - AM^s sin^ 9 + 8/3M/Visin 29 

~ l + cos2^ + 4MVssin2^ ' 

^ cos^ + 4/5M/v/isin^ 



'l + cos2^ + 4MVssin2^ ' 
Dnn = Dss = 0, 

.2. 1-4MVs 

Ann — Cnn = — sin 6* — -ttt'/ — ^^TT^ ■ 

1 + cos2 + AAP/s sm^ 

_ _ . sm0{l + 4Mys) - 16/3M/^cosg 
Css - -sm9 i + cos2^ + 4MVssin2^ ' 

_ _^M/^sm9cos9 - (3cos'^9 + A(3M'^/ssm'^9 
~ l + cos2^ + 4M7ssin2e ' 

Cln — Cns — . (87) 



As for single spin asymmetries (of. the remark in Sect. 3.1) the observables Dnn, Cns 
and Cln are related to the imaginary parts of interference terms and are therefore zero 
to lowest order of perturbative QCD. To higher orders they may be non-zero. In Fig. 7 
we show results for some of the spin correlations. Since these predictions extend to c.m.s. 
scattering angles up to 90° the /3 dependence of the spin correlations has been fully taken 
into account (beyond linear order) in the numerical evaluation. 



6 Summary 

We have analyzed the exclusive production of A+A~ pairs in proton-antiproton collisions. 
We have argued that this process is mediated by collinear emission of soft u (u) quarks 
from the (anti-)proton, a hard scattering uu —>■ cc and subsequently a collinear absorption 

*0f course, the subprocess amplitudes have to be transformed to the c.m.s. basis analogously to (57) 
and (58). 
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Figure 7: Final state helicity correlations Cij and the p — Ac polarization transfer Dn 
versus cos^ at s = 30 GeV^. (colors online) 

of c (c) quarks by the remainders of the (anti-)proton turning them into the final 
(A^). The hard subprocess is calculable within perturbative QCD to a given order. As 
we showed the hadronic matrix elements that describe the p —>■ A^ {p —>■ A^) transition, 
can be parameterized in terms of GPDs. The process amplitudes thus factorize into a 
product of hard subprocess amplitudes and 1/x moments of GPDs, i.e. generalized form 
factors as, for instance, in wide-angle Compton scattering. 

To model the GPDs we employ the representation of the GPDs as overlaps of LCWFs 
for the involved baryons. The LCWFs for the valence Fock states of the baryons are 
taken from the literature. Although, in general, the restriction to the valence Fock states 
is insufficient for modelling the GPDs as a function of three variables x, ^ and t, it is 
expected to be a reasonable approximation for our case. As is predicted by the HQET, 
the Ac valence Fock state is dominated by the simple configuration where the c quark 
carries the helicity of the Ac while the light quarks are coupled in a spin and isospin zero 
state with only little admixture of other configurations. It is also expected that higher 
Fock states play only a minor role for the Ac. A possible more complicated structure of the 
proton involving, for instance, parton orbital angular momentum or higher Fock states, is 
irrelevant since its overlap with the simple Ac LCWF is zero. We therefore achieve rather 
robust model GPDs with only a few parameters, like the transverse size parameters or the 
normalizations of the LCWFs which can be adjusted within certain ranges without loosing 
their physical interpretation. With the GPDs at hand we are in the position to evaluate 
the form factors and, subsequently, to predict differential and integrated cross section as 
well as a number of spin correlation parameters. We emphasize that the prediction of the 
integrated cross section achieved here does not differ much from a previous estimate in a 
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model that bears resemblance with the handbag approach [5]. 

We also note that the u ^ c transition GPDs do not only form the soft physics input 
for the handbag contribution to the process pp — > A+A~ but allows also to calculate 
the — > p transition form factors. Another possible application is elcctroproduction of 
charmed mesons (e.g. ep e'D(Z)*)A+). Extensions to other charmed baryons and to 
calculation of corrections to the subprocess is left to a forthcoming paper. 
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Appendices 



A Kinematics 

The momenta and light-cone helicities of the incoming proton and antiproton are denoted 
by p, 11 and g, i/, those of the outgoing A+ and A~ by p\ /i', and q', u', respectively. The 
mass of the proton is denoted by m that of the A+ by M. We work in a center of mass 
frame where the baryon's momenta in light-cone coordinates are parameterized as 



p = 



Q 



2(1 + 0^' 



+ A5_/4 



p = 



M2 + Ai/4 
2p+(l - ' 



M2 + Al/4 

2r(i-o ' 



Op' 



2 
A, 



2(1 + 0p+' ^' ■ ' 2 
It is convenient to introduce sum and differences of the baryon momenta 

P = 2p = p + p' , Q = Q -\- q' , A = p'— p = q — q' 

The momentum transfer A is expressed by 



A = 



(88) 



(89) 



(90) 



2p+(i-a 

The three-components of the c.m.s. momenta of the incoming (anti) proton, (— )p3, and 
the outgoing (anti) Ac, (— )P3, read 



P3 



1 



(91) 



where s = (p + q)"^ — ip' + q'y is the usual Mandelstam energy variable. We also introduce 
the abbreviation 

Ai = Jl-Amj/s (92) 
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in a more generic notation (m^ stands for m, M, or the charm-quark mass mc, for rrii = rric 
Mandelstam s has to be replaced by s). While the three-component of the incoming 
proton's momentum is always positive that of the outgoing Ac can become negative for 
large c.m.s. scattering angles because of the unequal-mass kinematics. This change of 
signs occurs at 

Thus, due to the unequal-mass kinematics, is zero for forward scattering, reaches 
the maximal value and decreases again towards zero for backward scattering. Skewness, 
defined as a ratio of light-cone plus-components of baryon momenta, ^ = (p+ — ) / (p+ -|- 
p'"^), can be expressed as 



2 + ^Kl - Ai/s + sign(p'3)VAM-AiA 



A? 



A2 



.2 ^^2 



2A™Amcos6' 2 + ^/kl 



2s 



3 + + Al/2 
^2 s 



A2 



for s 



2Kmk.M COS 9 



oo 



(94) 



where 9 is the c.m.s. scattering angle. In a similar fashion can the other parameter 
occurring in the definition of the baryon momenta (88) be written as 



A^-Ai/. + sign(p;,)y'Al,-Ai/. 
^y|[2 + ^A^ + A2^ + 2A^AMCOse] . 



(95) 



As a consequence of the unequal mass kinematics ^ cannot become zero. For p'^ > 0, 
however, the skewness is fairly small and tends to zero for s ^ oo. 
The squared invariant momentum transfer is given by 



t 



A^ = 
A2 



2e 



1-e 

s 



1-e 



A^ + Al - 2sign(p'3)^A^ - Al/s^A^ - Ai/.] 
- 2A^Am cos 9 



It cannot become zero for forward scattering but acquires the value (^o = ^( A 
0)): 



± 



(96) 



to = i(Ai = 0,p^>0) 



2^0 



(l + 6)M2-(l-eo)m^ 



--(A^-Am)'. 



(97) 
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It is convenient to introduce a variable t' that vanishes for forward scattering. It is defined 

by 



t' 



S 

1 — COS 9 



A^Am - sign(p^)^A^ - ^l/s^Kl, - Ai/s] 



(98) 



Solving (98) for A^, one finds 

sA„iAm + 1' 



A^ = -t' 



s/4(A^ + KmY + t' A'i + Al^ + 2A^Am cos 9 ' 



(99) 



Sine and cosine of half the c.m.s. scattering angle 9 are connected with the Mandelstam 
variables by 



sin2^/2 = 
cos2^/2 = 



1 — cos 9 to — t 



2 sA^Am ' 

1 + COS^ Ui — u 

2 sA^Am ' 



where 



u, = u{Al = 0,p^ < 0) = -- (A^ - Am)' 



For forward scattering u reads 



uo = u{Al = 0,p'3 > 0) = -- (A„ + Am)' 



B The u ^ c generalized parton distributions 

The vector current of bilocal quark field operators is defined as 

V''{-z-/2,z-/2) = ^'^(-z-/2)7'^*"(^-/2) - *"(2-/2)7'**'^(-^-/2) . 
The Fourier transform of the plus component of its transition matrix element 



n-/cu _ ^ 



27r 



.e^^^""^-(A+;pV I V+{-z-/2,z-/2) \ p;ppi) , 



can be decomposed into five covariant structures 



M + m 



u{p')u{p) , 



A+ 



M + m 



u{p')u{p) , u{p')-f+u{p) , 



M + m 



M + m 



(100) 
(101) 

(102) 



(103) 



(104) 



(105) 
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The combinations of momenta are defined in (89). Helicity labels are omitted wherever it 
can be done without loss legibihty. Only two of the covariant structures are independent 
since we have the familiar Gordon decomposition and analogous relations at disposal 



Pf'u{p')u{p) = (M + m)u{p'Yi^'u{p) - m{p')af"'A,,u{p) 
A^m(p')m(p) = (M - m)u{p')j^u{p) - tu{p')a'"' Puu{p) . 

In addition we have for the light-cone projections the relation 



(106) 



(107) 



see (88) and (89). The two independent GPDs are chosen in analogy to the flavor-diagonal 
case (cf. [48]) 



cu 



u{p',iJ,') 



M + m 



u{p,ij) 



(108) 



In the unequal mass-case the matrix elements of the local vector current describing 
the electroweak p A+ transition can be decomposed into three covariant structures for 
which one may choose 



(A+;py I ^^^^(0,0) I = iZ(p',//') Fy^^^ + Ft 



m + M 



+ F. 



M + mi 



uip,fj,). (109) 



The four- vectors and P'* are independent; only their plus components are proportional 
to each other. 

Let us integrate (104) and (108) with respect to x, which reduces the bilocal matrix 
element to a local one that describes the weak transition form factors for which we insert 
the form factor decomposition (109) 



dx'H.^i^ 



(A+;pV| V+(0,0)|p;p/x) 



= y ^ dxH'^^x, t) u{p\ i^'h+u{p, dxE'^^x, t) u{p\ i^')^^^u{p, /.) 



u{p', /i') 



?rr+''A A 



m + M 



M + m 



M + m 
u{p,ii). (110) 



Using again the Gordon decomposition (106) for the term oc F3, we find the following 
sum rules 



J ^dxH^^{x,U) = Fv{t) - CFsit) , 
j ME''^{x,i,t) = FT{t) + iF,{t) . 



(Ill) 
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In the limit of equal hadron masses conservation of the vector current requires F3 = 0. To 
lowest order of the electroweak theory the weak transition form factors are zero; this holds 
in particular at the scale of Mz- For lower scales, like M, there are effective operators 
that mediate u ^ c transitions through finite Wilson coefficients. 

The lower vertex is treated analogously. It is related to the upper one by charge 
conjugation. It is then easy to show that 



(A+ I V-^i-z/2,z/2) I p) = -(A; I V+{-z/2,z/2) \ p) 



(112) 



Taking the Fourier transform of this matrix element in analogy to (104), making the 
replacements 

p^q, p'^q', A^-A, T^^^T", (113) 

and rewriting the right-hand side of (108) with the help of the behavior of the Dirac 
spinors under charge conjugation one finds for the lower vertex 



UCU - 

= v{q,u 



27r 



^rxq ^+(A;;gV I V-{-z+/2,z+/2) \p;qu) 



M + m 

The lowest moment of 7if/„ leads to the vector-current matrix element 



v{q',u') 



(114) 



(A-;gV| V-(0,0)|p;gi/) 



-1 



dxH^'^ix, C, t) v{q, y)l~v{q\ z/') - / dxE''''{x, ^, t) v{q, u 



M + m 



v{q',u' 



v{q,iy) 



Using 



m+M ''M+m 



viq',u'). 



and the Gordon decomposition 

Q^v{q)v{q') = -{M + m)v{q)^^'v{q') -iv{q)a^''A^v{q') , 

we find exactly the same sum rules as in (111). 

Next we proceed analogously for the axial vector transition matrix element: 

A^'{-z-/2,z-/2) = *'^(-z-/2)7'^75*"(z-/2) - *"(^-/2)7'^75*^(-z-/2) , 

dz~ 



(115) 
(116) 
(117) 



'IJCU _ 

i^ix'n — P 



27r 



e^^^""^-(A+;pV I A+{-z-/2,z-/2) \ p-p^i) , 



(118) 
(119) 
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can again be decomposed into five covariant structures 



P+ 



M + m 



A+ 



M + m 



M + m 

We also have the two relations 



M + m 



(120) 



M(p')^(T'^"A,75M(p) = (M - m)ix(p')7^75M(p) - P+m(p')75«(p) , (121) 

and the projection (107). The two independent GPDs are chosen in analogy to the flavor- 
diagonal case 

^-(x,e,t)7+ + ^'"(^,e,^) 



Wc define the general local axial-vector current matrix element as 



l5u{p, n) 



(122) 



r A*" 
(A+;pV I ^^^(0,0) I p■p^x) = u{p\^x') Gai' + Gp—^ + G, 



m + M 



M + m 



lbu{p, n) . 
(123) 

This may be compared with a decomposition which is useful for applications of the HQET 
[16] . Integration of the GPD with respect to x leads to 



dxHi;^ = (A+;pV|>1'^(0,0) 



M^"(x, t) u{p\ pl')^-^^5u{p, pi) + / dxE'^'^ix, t) 



= u{p', n') 



GAit)i^ +Gp{t)-^ + Gs{t) 



M + m 

Using (107) we finally obtain the sum rules 



M + m 



M + m 

15U{P, IJ) . 



u{p', l^')lbu{p, n) 
(124) 



dxH'^{x,^,t) = GA{t), 



dxE^-{x,tt) = Gp{t)--Gs{t). 

Applying charge conjugation 

(A+ I ^+(-^/2,^/2) I p) = -{A- I ^+(-^/2,^/2) I p) , 
one obtains the parameterization of the lower vertex 

dz+ 



(125) 



(126) 



nj cu 

rf'ij'i. 



27r 



^^xq --(A-;^V I A-(-^+/2,^+/2 I p;qiy) 



= v{q,u) 



H 7 75 - £/ 



A-75 " 

M + m 



v{q',iy'). 



(127) 
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The lowest moment becomes 
-1 



dxW^ = (A;; I A+(0, 0) I p- qv) 

^+ J dxH'^v{q, i^h'^^viq', v') - J dxE'^'^viq, 

Q- 



v{q,u) 



GA{t)i- - Gp{t) 



M + m 



+ Gs{t) 



M + m 



N ^ 75 / , ,x 

J^)jT^v{q ,u) 

M+m 

Ibviq', v') . 



(128) 



One reads again off the sum rules (125). 

In an analogous way the helicity-ffip GPDs can be written as (for definitions see Ref. 
[23]): 



dz 



X (A+;pV I *'(-;2"/2)«a+%"(;2-/2) - ^"(^-/2)z(t+%'=(-^-/2) | p;pi^) 



= u{p',ij!) 



+ E^l ^- + E^ 



2mM 

^+ pj — p+y 



M+m ^ M+m 

Charge conjugation symmetry leads to 

dz+ 



u{p,ij) 



(129) 



/ 



27r 



X (A7; I *^(-^+/2)z(7-^*"(z+/2) - *"(z+/2)z(7-^*'^(-z+/2) | p- qv) 



v{q,iy) 



cu 
T 



2mM 



+ 



M + m 



771 cu 
Jjjrp 



M + m 



v{q',iy'). 



With the relations 



a^Ss = +^a^' 



(130) 



(131) 



a combination of cr"*"^ and cr"*"^ matrix elements (129) and (130) can be reexpressed as 

(132) 
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